Under consideration for publication zn J. Fluid Meek. 


1 


Very, very fast wetting 


By DAVID JACQMIN 

NASA Glenn Research Center, Cleveland, OH 44135, USA 


(Received 30 November 2001) 

Just after formation, optical fibers aie wetted stably with acrylate at capillary numbers 
routinely exceeding 1000. It is hypothesized that this is possible because of dissolution 
of air into the liquid coating. A lubrication/boundary integral analysis that includes gas 
diffusion and solubility is developed. It is applied using conservatively estimated solubility 
and diffusivity coefficients and solutions are found that are consistent with industry 
practice and with the hypothesis. The results also agree with the claim of Deneka, Kar 
k Mensah (1988) that the use of high solubility gases to bathe a wetting line allows 
significantly greater wetting speeds. The solutions indicate a maximum speed of wetting 
which increases with gas solubility and with reduction in wetting-channel diameter. 
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1. Introduction 

Immediately after formation, optical fibers, in order to improve their resiliency and 
refractive properties, are wetted and then coated by running through a polymer bath. 
This wetting typically takes place at capillary numbers Ca = of 0(1000), that is, 

at wetting speeds U 0 of about 1000 to 2000 centimeters a second, surface tensions a of 
about 20 to 30 dynes/cm, and liquid viscosities pi of about 20 to 40 poise (Dimitropoulos 
et al., 2000; Ravinutalaet al., 2000; Jochem k Ligt, 1985; Lyytikainen, 1998). 

Optical fiber wetting is somewhat like plunge-tank wetting; in both the material being 
wetted is drawn through a surrounding free surface. However, with optical fibers wetting 
takes place in a very narrow orifice - the entrance channel or die of the coating apparatus 
- that is typically only 300 to 600 microns in diameter (Jochem k Ligt; Dimitropoulos 
et al.; Ravinutala et al.). The fiber itself has a diameter of about 125 microns. The gap 
between the fiber and the orifice wall is thus typically about 100 to 300 microns. The 
drag of the fiber on the liquid is resisted by the application of high pressures that force- 
feed the liquid up into the orifice. Because of the fiber’s high speed and the relatively 
high viscosity of the fluid, pressure gradients in the entrance channel are on the order 
of atmospheres per millimeter. It appears that proper shaping of the inlet directs these 
pressures to the wetting line so that air entrainment can be resisted. 

These high pressure gradients are undoubtedly helpful in resisting air entrainment. 
The small size and axisymmetry of the fiber may also be helpful. Simpkins k Kuck 
(2000) reported a critical capillary number (the maximum Ca for successful resistance 
of air entrainment) of 2.1 for small-diameter fibers entering unpressurized glycerin. By 
comparison, for flat tapes in plunge-tank experiments critical Ca is less than 1. Tor 
pressurized orifices Simpkins k Kuck have been able to observe and photograph successful 
wetting at Ca greater than 20 (private communication). 

However, the success of high speed optical fiber wetting cannot be entirely explained in 
terms of global quantities. Because of the high wetting speed the liquid dynamic contact 
angle must be at or near 180°. As will be discussed below, the potential exists for very 
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high and destabilizing pressures to be built up in the receding gas phase near the wetting 
line. Some theory is needed for how these small-scale but very large stresses are relieved 
or resolved. 

A hint of how this is accomplished is given by the frequent practice in the industry of 
bathing the wetting line with high solubility gases. The efficacy of this was first claimed 
in a patent assigned to Corning Glassworks (Deneka, Kar & Mensah, 1988). The gas 
recommended and claimed in the Corning patent was carbon dioxide, which is typically 
3 to 5 times as soluble in polymers as is air. Jochem & Ligt (1985) experimented with 
CCI 2 F 2 , which is even more soluble in polymers than CO 2 , and showed great improvement 
in coating speed and quality. With CCI 2 F 2 they demonstrated bubble-free coating at 
1250 centimeters per second. With air, using the same wetting die, bubble entrainment 
occurred at 300 cm/sec (Jochem Sz Ligt, 1987). 

The hypothesis indicated from these claims and results is that dissolution of gas into 
t he wetting fluid serves to significantly relieve high gas pressures, thereby facili tatin g high 
$peed wettin g. I n orderto examine this hypothesis, this paper disc usses a combination 
lubrication /boundary integral analysis that quantifies gas solubility effects. The analysis 
approximates the dynamic contact angle as being 180°, giving a gas column shape that 
forms a half cusp between the wall and the liquid. 

Gas properties that could have important effects on high-speed wetting stability include 
1) viscosity, 2) compressibility, 3) solubility in the liquid, 4) diffusivity in the liquid, 5) slip 
along the fiber, 6) Knudsen diffusion of momentum, significant for gas column thickness 
less than about 300 nanometers, and 7) increased local gas solubility when the gas phase 
is confined to a very small domain (nanoscale). Properties of the liquid phase include 1) 
viscosity, 2) surface tension, 3) van der Waals attractions between the liquid and the solid 
and 4) decrease in the surface tension as the liquid aproaches the solid. The main body 
of this paper considers the first four gas properties and the first two of the liquid. An 
appendix briefly discusses the modeling and effects of most of the remaining properties. 
Also, Jacqmin (2001) gives a more complete description of modeling of the gas phase. 
All the properties relegated to the appendix are stabilizing. Their neglect in the main 
article thus permits a relatively streamlined analysis of the “worst case” . 

In spite of the huge practical and fundamental scientific interest of high-speed wetting, 
theoretical attempts aimed at understanding it are largely lacking. An early significant 
result is by Benney & Timson (1980), who derived an eigenfunction solution for 180° 
wetting of a solid. Their solution gives the height H(x) of the liquid-gas interface (the 
distance of the interface from the fiber or, equivalently, the gas-column width) as behaving 
like |x| m+1 , where m, the eigenvalue, equals — arctan(2Ca)/7r.f As Ca goes to oo, m 
approaches 1/2 from above. The solution neglects the effect of gas stresses. 

There has been more work on free-surface cusp flows (which can be viewed as “roll-on” 
liquid-liquid wetting). Joseph et al. (1991), showed that these have the same eigenfunc- 
tion solution as the one found by Benney fc Timson. Shortly after, Jeong & Moffat (1992) 
derived exact solutions for near-cusp flows using complex-analysis techniques. These so- 
lutions also neglect gas and wetting stresses. At micro- and mesoscale distances from 
the cusp these solutions approximately match to the Benney fc Timson eigenfunction. 
Jeong & Moffat showed that the coefficient of the eigenfunction is determined by the 
outer length scale. In their case this scale was the depth of a vortex dipole that drove 
their flow. At the nanoscale their solutions show a slight rounding of the cusp, which 
eliminates the singularity in surface tension forcing that would otherwise seem to occur 
(a delta function with amplitude 2<r). 

f This corrects a sign error, introduced at their equation (2.12) and carrying on to (3.9). 


3 


Very, very fast wetting 

The difficulty with either cusp or roll-on flow is that the pressure in the displaced 
gas phase can become very high, leading to instabilities and gas entrainment. The gas 
phase is swept into the cusp region by the movement of the fluid and of the solid being 
wetted. Since most of the gas has to return, a Poiseuill e-like flow is established with 
pressure gradients like 1 / H 2 . If the flow were precisely cusp or roll-on, and if the gas 
were completely insoluble in the fluid and obeyed the Navier-Stokes equations, then the 
gas pressure would be fiercely singular; for the Benney & Timson solution it would be 
like This would dwarf the viscous and surface-tension-related stresses in the 

liquid. Shikhmurzaev (1998) speculated that this problem might be relieved by Knudsen 
diffusion. However, this also leads to a singularity (like x~ m , see Jacqmin). Nor does slip 
along the wall solve the difficulty, unless slip is also allowed along the liquid-gas interface. 
Eggers (2001) recently conducted a lubrication/boundary integral analysis of insoluble 
two-phase cusp flow assuming, following Jeong & Moffat, that the cusp is actually slightly 
rounded. He estimated that failure occurs for Xk 3 ^ > r, where A is the viscosity ratio 
/i gas //ii iq , k is the curvature at the rounded cusp, and r is a to-be-determined 0(1) 
quantity. Since, from Jeong & Moffat, k — C\ exp(2 7rCa)+ where c\ is 0(1), the maximum 
capillary number for cusp stability would be 

Ca maK ~ A h 2- + 0(1) (1.1) 

O 7T A 

This is way too low to account for the stability of optical fiber wetting. 

Other industrial processes besides optical fiber manufacturing are reported to allow 
high Ca coating. For example, curtain coating (where a metered stream of liquid falls 
from a height above a moving substrate and then coats the substrate) allows wetting up 
to a Ca of about 10 (Kistler, 1993, p. 343). It seems certain that the additional force and 
pressures associated with the fall of the liquid act to allow the observed higher wetting 
speed. Recent experiments by Blake, Bracke & Shikhmurzaev (1999) have quantified 
some aspects of these effects. Kistler (1993, pp. 339-346) gives a thorough review of the 
various possible macro- and micro-scale phenomena that may affect air entrainment and 
wetting stability. 

The following first derives, for given H{x ), an integro- differential equation for the 
pressure p g in the gas phase. This uses a lubrication approximation for the gas flow and 
a Green’s function analysis for the diffusion of the gas into the liquid. H{x) is then found 
in terms of an integral of the stresses on the liquid interface. This part of the analysis is 
linearized about a flat interface. Since the interfacial stresses can be expressed in terms of 
H and p & the system of two equations is closed. We find that there are both minimum and 
maximum speeds for which successful wetting is possible. For high capillary number these 
cutoff speeds are expressible in terms of a function of two nondimensional parameters. 
The equations are applied to the regime typical of optical fiber coating. It is shown that, 
indeed, gas solubility significantly mitigates gas pressure and that, in agreement with 
the claim of Deneka et al., the use of higher-solubility gases makes high wetting speeds 
easier to achieve. 


2. The gas flow 

We consider the flow of an isothermal, soluble, ideal gas in a “half-cusp” region bounded 
by a solid wall located at y = 0 and a liquid interface at y — H{x). The half cusp is in 
x > 0 and ends at the wetting line at x — 0, see figure 1. The bounding surfaces are 

% This uses the definition of Ca employed by Joseph et al., which is the same as used in this 
paper and 16 times that of Jeong &; Moffat’s. 
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Figure 1. The wetting geometry, x is the horizontal coordinate (increasing toward the 
right), y the vertical. The wetting line is at x — 0, y = 0. The gas is thus confined to 
x > 0. The solid is moving to the left with velocity — Co. 


quasi-parallel, allowing a lubrication flow analysis for the gas. The solid is moving at 
velocity — Co and the liquid motion is approximated as also being — Co* The lubrication 
equations for the gas flow are 


duPg fcPg _ 0 
dx dy 


dpg _ & 2u 

dx d y 2 


Pg = PgO) 


(2.1) 


where u and v are the gas x- and y - direction velocity components, pg is the gas pressure, 
and is the gas viscosity. The first equation is continuity, the second, z-momentum. 
The convection-diffusion equation for the gas in the liquid can be written as 


■V.g.D\ 

OX 


d 2 Q d 2 c { 

dx 2 + dy 2 


(2.2) 


Ci is the gas concentration in the liquid and D is the diffusivity. Ci is assigned units of 
pressure. 

The boundary conditions at the liquid-gas interface are 1) no-slip, u = -Co, 2) phase 
equilibrium, = Spg(H~) } where S is the Henry’s law solubility coefficient for 


the gas in the liquid, and 3) equality of normal fluxes, —D 


dy 


h+ 


— (Ug as ^liq)Pg 




Boundary conditions at the gas-wall interface are 1) no normal velocity and 2) no-slip. 
The gas is at atmospheric pressure, designated p OQ) at x = Too. For simplicity, it will be 
assumed that the fluid is saturated with the gas {C\ equals S times Poo) at x — Too and 
y = Too. 

Integration of the continuity equation from 0 to H{x) gives 


dU^Hpg 

dx 


— (^gas ^interface )Pg 


— ^diffusion Pg 


H~ 


= —D 


H- 


dC i 
dy 


H+ 


[7 av = —H~ l f udy. From the momentum equation, 

1 H 2 dp R 

U 0 — U av — — — 77 r 


p g 12 dx 


(2.3) 


(2.4) 


Substituting (2.4) into (2.3) yields 

ff 3 dpg 
12p s Pe dx 


+ fa\ U ° Hp z 




m 

dy 


H+ 


(2.5) 
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In the convection-diffusion equation z-transport is dominated by convection. The in- 
direction diffusion term can be neglected and the equation reduces to 

d 2 Cj 


- U 0 = D- 
ox oy z 


(2.6) 


dCi/dy also obeys this equation. A Green’s function solution is available (Carslaw & 
Jaeger, 1959f) that expresses dC\/dy in terms of the boundary values of d 2 C\/ dy 2 : 


dOi 

dy 


D 

ttUo 


f 


d 2 Q 


dy 2 


dx' 

e 4D(x'-*) ■ 

= /f+ v x * - 


This, from (2.6), can be rewritten as 



_ Un(y~H) J 
g 4Z)(x'-r) 


dx' 


y=H+ 


y/x 


(2.7) 


( 2 . 8 ) 


Equation (2.8) can be used to find the flux at the interface in terms of the cusp gas 

dCi \ 


pressure. Multiplying by D , using 
to y — H + gives 


D 


dQ 

dy 


= S 


//+ 


y=H+ 


DU 0 


- Q^Es 

b dx' 


and specializing the equation 


y=H- 


r dp, 

Jx dx' 


dx' 


H- 


y/x' — X 


Applying this to equation (2.5) then gives 


dx \ I2p, 


R3 v 

g dx 


d- 


±(UoHp.) = l 


00 dp, dx' 
dx' yj x* — x 


(2.9) 


(2.10) 


a one-dimensional integro-differential equation for the steady-state pressure of the gas in 
the half cusp. 


T 

3. The liquid flow 

The flow in the liquid is taken to be Stokes flow. (jTo leading orde£ its velocity is — Uq. 
Linearizing about — Uo and a flat interface, the ki nematic equ atlon for the free surface 
is, 

-Uo^- = vj (3.1) 

OX 

where vj is the normal liquid velocity at the interface, vj is related to the interfacial 
stresses by the boundary integral 


1 

27 r/q 



x\ fydx' 


(3.2) 


(Pozrikides, 1992) where f y is to first order the interfacial normal stress minus the at- 
mospheric pressure. To this order shear stress contributions to vj are negligible. Taking 
the first derivative of (3.2) gives 


= - J- r -A*’ 

dx 2? r/jt\ J_ (X) x 1 - x 


(3.3) 


f See equation (9) in §2,9. The t occuring in the argument of the exponential function should 
be r. 
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a Hilbert transform of f y j 2p\. It can be inverted to 

f ~ r — — 

V n J-oo dx' x 1 — x 

and, since vj in x < 0 is 0, specialized to 

2 fi\ f°° dvj dx' 

V 7T 0x f X f — X 


(3.4) 


(3.5) 


f y — <jH xx - Poo and dvj/dx = — U$H XX . Equation (3.5) can thus be arranged in 
the form of an integral equation for H xx forced by p g : 

nr 2/Wo f°° H X ' X . J , , Q ^ 

<rH xx / — dx = Pec - P g (3.0) 

7 T J Q X f -X 

The equation is linear with a singular Cauchy kernel. The solution for H XX} as discussed 
by Mikhlin (1957, pp. 126-131), is 


d 2 H 

dx 2 


(pg ~ Poo)/* 
1 + 4Ca 2 


2 Ca/rtr m 
1 + 4 Ca 2 



(z') 1 ~ m 
x f — X 


{Pz{ x ')-Poo)dx' + c 00 x m 1 


(3.7) 


Equation (3.6) corresponds to Mikhlin’s equation (1), p. 127, (3.7) corresponds to (23). 
Coo is a free parameter. C 00 x rn ^ l i where m — — arctan(2Ca)/7r, is the homogenous 
solution to the unforced problem. It is the eigenfunction first found by Benney & Timson. 
The arctan branch is taken such that m varies from 1 (as Ca — ► 0) to 1/2 (as Ca — ► oo). 
This is the only branch that, in equation (3.6), gives a finite integral. 

The advantage of equation (3.7) is that p g - p^ is significant only in the microscale 
region near the wetting line whereas H xx in (3.6) decays very slowly. As x — ■ * oo the 
first term on the right hand side of (3.7) decays like a: -Ki+ 2m ) j the second like x m ~ 2 and 
the eigenfunction term like x m ~ l . Thus the Benney & Timson eigenfunction becomes 
dominant and gives the macroscopic shape of the fluid interface. Coo is therefore the pa- 
rameter that allows connection of the macroscale appearance of the flow to its microscale 
behavior. Ultimately, as discussed by Jeong & Moffat (p. 11, their variable c), its value 
is related to the outer length scale set by the flow geometry. Since p g is nonsingular, H xx 
also behaves like x m ~ 1 as x — ► 0. The coefficient of x m ~ l there is 


c 0 = jTVr m M*') - 


(3.8) 


4. The dimensional equations. Ranges of parameters 

The coupled system of integro-differential equations is thus: 


d ( H 3 dp d / _ a _ /DUo dp g dx' , A 

dx\ 12/z g Ps dx ) + dx \ UoHp * ) " 1 r I dx 1 ^ ^ 


(*o 


t\l — m 


-(Pg( x ')-Poc,)dx'+CooX m 1 (4.2) 


d 2 H _ (Pg — Poo)/* 2Ca/ir(T 
Ix 2 ^ 1 + 4Ca 2 1 + 4Ca 2: Jo 

The boundary conditions for the equations are that p g is finite and H = H x = 0 at x = 0 
and that p g equals at x — oo. The far-field behavior of p g is 


-■/ 


dpg ^ 12pgUo 12p g 

dx H 2 i7 3 p g 


(4.3) 


7 


Very , very fast wetting 

where L is the total rate of dissolution of air into the fluid, the integral from 0 to oo of 
the right hand side term in (4.1). It is found as part of the solution. The integral from 0 
to oo of the first two terms of equation (4.2) is precisely 0. (This result can be shown by 
converting the integral in x in the second term to a Mellin transform in x/x'.) Because 
of this the leading order behavior of H at oo is of the form CooX m+l /m(m + 1) pl us a 
constant times x m . 

The equations have seven independent parameters, Co, /i g} Poo, SVD, a, fi\ (in Ca ), and 
Coo We are primarily interested in the effects of varying S, Coo, and Co- In optical fiber 
coating the other parameters show less significant variation. A single set of representative 
values of p g , /q, D y <r> and has therefore been used for all dimensional calculations. 
In cgs units, this set is p g — .00019 poise, /q = 30 poise, D = .00003 cm 2 /sec, a — 30 
dynes/cm, and poo = 1.01 x 10 6 dynes/cm 2 (one atmosphere). Durrill k Griskey (1966, 
1969) give D for a number of different gases diffusing in various molten polymers. The 
value chosen for these calculations is in the middle range of their data. 

Durrill k Griskey also give gas solubilities. From them, carbon dioxide is typically 
about 4 times as soluble in molten polymers as is air. A typical value of 5 for air is about 
.1, for CO 2 about .4.| Both cases will be considered. 

Following the results of Jeong k Moffat, it will be assumed that Coo i s approximately 
equal to the inverse square- root of the lengthscale of the outer geometry. In their case 
the outer length scale is vortex dipole depth d and they find a Coo of 1.225 fy/d. The 
relevant outer length scale for optical fiber wetting is the gap between the optical fiber 
and the wetting channel’s wall. This is typically 100-300 microns, indicating a range for 
Coo °f about 5 to 10. 

Optical fiber coating speeds mentioned in the patent literature range from about 300 
cm/sec to over 2000. Corresponding capillary numbers range from over 100 to over 2000. 
The industry norm is in about the middle of this range. These very high capillary numbers 
will allow some useful simplification of the analysis. 


5. Existence of solutions 

To be physically valid, solutions for H must be positive, but p g , which is also necessarily 
positive, tends to make H xx and thus H negative. At the wetting line, negative H xx would 
immediately cause H to be negative. Co must therefore be greater than 0 and so, from 
(3.8), it must be true that 

a icii r^~ m (5-i) 

For large Ca a cutoff value of Coo exists below which there is no solution. This value is 
a function of the seven equation parameters. 

At large Ca, for given /i g , pco, SVD } a, /q, and Coo, there can be either no Co that 
allows solutions or a range of solutions from a minimum Co (a lower Co cutoff) to a 
maximum (the upper cutoff). Solutions also exist as Ca — ► 0 in a small region near and 
below Ca — 1. 

A simple class of solutions exists for large enough Coo. If H were of the form otx^ then 
the shape of p g and its maximum value would be independent of a. This is because a can 
then be transformed out of (4.1) by the substitution x — ► At large Coo H 

approximates this form because the third term on the right hand side of (4.2) dominates 

f This paper uses a somewhat different definition of the coefficient, gas density in the polymer 
divided by the density of the gas, than do Durrill k Griskey. This increases the coefficients given 
by Durrill k Griskey by a factor of about 1.6. 
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the first two. As Coo — * oo the first two terms on the right hand side of (4.2) stay finite 
and become independent of Coo- (4.1) and (4.2) thus partially decouple. 

The lower Uq cutoff occurs where p g — poo is small but at high enough velocity so that 
the coefficient of the integral in (4.2) is still O(l/{7o). At small p g — poo (4.1) can be 
simplified to 


d ( # 3 dp^\ j , H \ = _S_ I DU o f°° dp, dx' 

dx\ 12/Zg dx J dx \ 0 / Poo V 7T J x dx f yjx 1 — x 


(5.2) 


Scaling H like CocL™* 1 , where L s is the inner length scale, and finding L s by balancing 
the first and third terms in (5.2) and then the pressure scale p s by balancing the second 
and third, one finds that p g - Poo scales like Uq /3 . The integral term in (4.2) scales like 

(p g — Poo) /Co and therefore like Uq 1 ^ 3 , thus indicating the cutoff. The essential point 
is that the cutoff occurs in a region where a diminishment in Co causes a less rapid 
diminishment in pressure forcing and so, see equat,ion(3.6) with the comparatively small 
(j H xx term neglected, the perturbation to H caused by the pressure forcing increases. 


6. Nondimensionalization. Specialization to Ca >> 1 

The seven parameters can be reduced to three through nondimensionalization. Setting 

x^L,x , H^C oa L’P +1 H , p 6 -p 5 p 6 (6.1) 


where 


£m+l/2 


Co 


Ps - 


_ PgUj 


S 2 D 


(4. 1-4.2) become 

d f H 3 _ dp g \ dr - \ _ 1 f°° dp g dx^ 

dx \ 12 dx J dx v ^ s / y/r dx' yjx' — x 


( 6 . 2 ) 


(6.3) 


<Ph_ 

dx 2 


r (pg Pco/^Ps} 

1 + 4Ca 2 


‘lYCaj'K 
1 + 4 Ca? X 






(Pe( x ')- ^-)dx+x m 1 (6.4) 


x — X 


The resulting nondimensional parameters are Ca, Poc/ps and T = p s / <jC 0 oL r ?~ l • The 
boundary condition for p g at x = -foo becomes that p g approaches Poo/p® • Since Coo has 
been transformed to 1 the problem of finding its cutoff value has been changed to finding 
the boundary of the region in (Ca, poo/p* , T) for which solutions exist. 


As Ca — ► oo, m — ► 1/2 and the lengthscale L s becomes 


The first term on 


the right hand side of (6.4) is 0(1/Ca) compared to the second and can be neglected. 
Equation (6.4) reduces to 


d 2 H 
dx 2 


2lT pi 


' Uo \ 3 / 4 __ 
v C 2 S 2 D J X 


1/2 





(6.5) 


while (6.3) remains the same. The number of nondimensional parameters in the equations 
has been reduced to two, (p g /27rpi)(Co/C^S 2 D) 3 / 4 and Pco/Ps = S 7 Dp 00 /p^U 2 . Find- 
ing the solution region is now a matter of finding the cutoff value of the first parameter 
as a function of the second. 
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| S 2 D Poo /p g US | It)" 5 | 1(T 4 | It)' 3 | 10" 2 | .1 | 1. | 10. | 100. | 

| ( l t g /2nm)(U 0 /Cl o S 2 D) 3/4 1 4.560 | 2.531 1 1.361 [ 0.687 | 0.314 1 0.129 | 0.048 | 0.016 | 


Table 1. Cutoff values of (p g /2nfi])(Uo/C^ 0 S 2 D ) 3 ^ 4 for given S 2 Dpoo / . 


7. Results 

The equations have been solved numerically using spatially varying node-point separa- 
tions. For dimensional calculations, grid spacing was typically set to nanometer lengths 
near the wetting line, up to micron lengths in the outer region. A similar variation in spac- 
ing was used for nondimensional calculations. The derivatives in (4.1) were discretized 
using central differences. The integral was solved by approximating dpg/dx' as piecewise 
constant. The result can be integrated analytically. The integral in (4.2) was approxi- 
mated by taking (p g (x')-p 0 0 )(x') 1 “ m as piecewise linear. The resulting integrand is then 
also analytically integrable. Convergence checks were carried out by varying the number 
of points, grid stretching parameters, and the length of the calculated region. The nu- 
merical solution at the outer calculated point was fitted to the analytic far-field solution. 
Equation (4.1) can display singular behavior at x ~ 0. This was avoided by solving it 
iteratively using a time-like approach. A dp^/ dt-Yike term was placed on the equation’s 
left hand side. The equation’s second derivative then acts like a diffusion operator and 
convergence to a finite solution is obtained without much difficulty. 

Table 1 gives calculated cutoff values of (fx g /27rp))(Uo/C% 0 S‘ 2 D) 3 / 4 as a function of 
S 2 Dp od /p s Uq . The relationship between the two is well approximated by s = -.89 - 
,41r - .023r 2 , where s and r are logs to the base 10 of the two variables. Figures 2 and 3 
give some dimensional results drawn from this approximation. They show cutoff values 
of Uq as a function of Coo for air and for CO 2 . The region of existence of solutions is 
between the upper and lower cutoffs. The values of pi, D, <r, and poo that are used 
aTe given in §4. The minimum value of Coo for CO 2 is 2.23, with Uq equaling 508 cm/sec. 
The minimum value of Coo for air is 4.47, with Uq of 116. 

Dimensonal calculations were also made. Of prime interest is maximum pressure de- 
veloped, to check on the assumptions of gas ideality and Henry’s Law, and of amount 
of gas absorbed by the liquid, to check on the possibility of later nucleation of bubbles. 
Calculations were made for Uq = 1000 cm/sec for Coo equal to 6, 8, and 10. Using air, the 
maximum gas pressure developed is, respectively, 147, 181, and 196 atmospheres. The 
region of high pressure is small; for Coo = 6, for example, p g > 100 atmospheres extends 
to only 75 nanometers from the wetting line, p g > 10 extends to about 300 nanometers, 
and p g ~Poo > 1 extends to about 1.2 microns. The maximum pressure decreases as Coo 
approaches cutoff, because of elongation of the air half-cusp and the resulting increase in 
liquid surface area that the gas can dissolve into. Consistent with this, the volume flux 
of air into the coating decreases with increasing Coo; it is 6.8 microliters a second per 
centimeter (in the direction perpendicular to x and y) for Coo — 6, 5.9 for Coo — 8) and 
5.3 for Coo = 10 (volumes at atmospheric pressure and 20° C). Since the fiber is moving 
at 1000 cm/sec this is equivalent to absorbing a layer of air only 53 to 68 nanometers 
thick. If there is later nucleation of bubbles they would be nanoscopic in scale. 

The same calculations for C0 2 yield 30.7, 31.8, and 32.3 atmospheres, with C0 2 fluxes 
into the coating of 7.74, 6.67, and 5.95 microliters per second-centimeter. The case Coo — 
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Figure 2. Cutoff values of Co as a function of Coo, for air. 5 - .1, = .00019 poise, 

- 30 poise, D = .00003 cm 2 /sec, a = 30 dynes/cm, and = 1.01 x 10 6 dynes/cm 2 . 



Figure 3. Cutoff values of Uq as a function of Coo, f° r CO 2 . S = .4, — .00019 poise, 

- 30 poise, D = .00003 cm 2 /sec, cr = 30 dynes/cm, and poo = 1.01 x 10 6 dynes/cm 2 . 

4 gives a maximum pressure of 28.2 and a flux of 9.58. The use of CO 2 instead of air 
greatly lessens the maximum gas pressure and allows a much smaller Coo* 

Figures 4 and 5 show pressure and gas solution rate as functions of x. The calculation 
is for air with U 0 = 1000 cm/sec and Coo = 5.3. The flow is close to the cutoff point. 
Of interest is that the pressure reaches a maximum away from the wetting line. The 
rate of dissolution is also a maximum there. Very near the wetting line the gas flow is 
determined by a balance between (4.1)*s second term and its integral. The gas in this 
region is in nearly plug flow. 


8. Conclusions 

This paper has attempted to explain the phenomenon of very high capillary number 
wetting that is common in the optical fiber coating industry. The usual model of wetting, 
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Figure 5. Gas solution rate as a function of x for air with Uq = 1000, Coo = 5.3. In 
milliliters at standard temperature and pressure per sec cm 2 . 


in which the two fluid phases are completely insoluble, points inevitably to wetting failure 
at an 0(1) capillary number (Eggers). However, standard practice in the optical fiber 
industry is to wet and coat fibers at Ca = 0(1000). It has been proposed here that 
stresses in the receding gas phase that would otherwise be destabilizing are significantly 
ameliorated by solution of the gas into the wetting liquid. The resulting model of wetting 
dynamics yields an adjustable macroscopic parameter that has a cutoff value below which 
there can be no solution. The calculated value of this cutoff is consistent with common 
practice in the sizing of optical-fiber-coating entrance dies. Maximum and minimum 
wetting speeds have been found as a function of gas and liquid properties and wetting- 
channel size. Results indicate that reduction in wetting-channel diameter allows higher 
wetting speeds. Results also agree with the claim (Denka et al.) that it can be very 
advantageous to bathe a wetting line with high-solubility gases. 




12 


D - Jacqmin 

I am very grateful to Dr. Enrique Rame of the National Center of Microgravity Re- 
search for his help and for many interesting and useful discussions. 


Appendix A. Slip, Knudsen diffusion, and van der Waals forces 

Slip and Knudsen diffusion effects are discussed in Jacqmin. With wall slip, equation 
(2.4) becomes 

1 ff 2 (\ H A l \ dp & 


U 0 - U w = — 


(Al) 


/i g 2 \2 ff + \ ' ff + \ 37 dx 

The wall slip length A is roughly equal to the mean free path. This is about 70 nanometers 
at standard temperature and pressure and is proportional to l/p g . There is no slip along 
the liquid (the two fluids intermingle) and so Uq — 27 a v continues to be proportional 
to H 2 dp g /dx. Flux due to Knudsen diffusion, which becomes important below about 
ff = 5A, is proportional to Hdpg/dx. Both effects produce a significant amelioration of 
stresses in the gas. The reduction in maximum pressure is over 50%. This results in turn 
in about a 50% reduction in gas absorption by the liquid. 

Van der Waals forces are normal to the interface and thus can be included in the model 
analysis. From Tsraelachvili, the van der Waals forces are equal to the rate of change with 
ff of the sum of the liquid and solid surface tensions (the excess free energies). They can 
be expressed in the form 


F= - 


A 


ff 


6tt (H + H o) 4 


d<Th dcr$ 
dff + dff 


(A. 2) 


For non-pol ar liquids Hq should be set to 1.65 angstroms (Israelachvili, 1991). When ff 
is greater than about one nanometer F can be set to the more common form — A/6nH 3 . 
A is the Hamaker constant. The actual change in surface tensions (as versus the first 
derivative of their change) becomes important only very close to the wetting line; for ff 
equal to one nanometer, ox is still about 97% of its value at oo. 

Calculations have been performed including both van der Waals forces and surface 
tension variation. The variation in surface tension turns out to be unimportant. Van der 
Waals attractions are stabilizing, partly in the same way as is large Coo, in that they 
cause an increase in ff x in the vicinity of the wetting line. They also, of course, resist 
the gas pressure. The importance of this is currently under study. Early results suggest 
that van der Waals forces may operate at too small a length scale to play a major role 
in setting either interface shape or the cutoff value of Coo- 
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